Influence of Capillary Condensation on the Near-Critical Solvation Force 
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We argue that in a fluid, or magnet, confined by adsorbing walls which favour liquid, or (+) phase, 
the solvation (Casimir) force in the vicinity of the critical point is strongly influenced by capillary 
condensation which occurs below the bulk critical temperature T c . At T slightly below and above 
T c , a small bulk field h < 0, which favours gas, or (— ) phase, leads to residual condensation and a 
solvation force which is much more attractive (at the same large wall separation) than that found 
exactly at the critical point. Our predictions are supported by results obtained from density- matrix 
renormalization-group calculations in a two-dimensional Ising strip subject to identical surface fields. 
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Finite-size contributions to the free energy of a fluid 
confined between two parallel walls, separated by a dis- 
tance L, give rise to a force per unit area between the 
walls or an excess pressure which is termed the solvation 
force fsoiv(L) (]]]. This force is, essentially, that which 
is measured in the surface force apparatus Theory 
predicts that at the critical point of a fluid the solvation 
force becomes long ranged as a result of critical fluctua- 
tions H , a phenomenon which is a direct analog of the 
well-known Casimir effect in electromagnetism The 
existence of the long-ranged critical Casimir force should 
be common to all systems characterized by fluctuating 
quantities with external constraints || . As yet there has 
been no direct, unambiguous experimental verification 
of the critical Casimir effect in fluids ||, although re- 
cent experiments do provide indirect evidence for its ex- 
istence ||. One of the difficulties is that the predicted 
leading power law decay of the Casimir force at bulk crit- 
icality, f so iv(L) ~ k B T c Ai 2 (d - l)L~ d , as L — > oo, is for 
(bulk) spatial dimension d — 3, of the same form as the 
solvation force arising from dispersion forces. Moreover 
the amplitude in many systems may be much smaller 
than the corresponding Hamaker constant (10]. The 
Casimir amplitude A 12 is a universal number, however, 
its value depends on the type of boundary conditions 
imposed on the two walls Q. The most relevant for 
experiments on pure fluids or for binary mixtures are 
symmetry breaking boundary conditions, i.e., when the 
confining walls exert surface fields on molecules in the 
fluid. Here we exploit the mapping between fluids and 
the Ising model and consider Ising spin systems subject 
to identical surface fields h\ = hi> 0. For these systems 
fsoiv{L) is expected to be attractive for all thermody- 
namic states. In d = 2, the Casimir amplitude is known 
exactly An = A ++ = — n/48 ||, whereas in d — 3 the 
most recent theoretical result is -0.428 M and the Monte 
Carlo estimate is -0.35 [Q. In this Letter we show that 
other, near-critical, thermodynamic states exhibit a sig- 



nificantly stronger solvation force (for fixed, large L) than 
that found exactly at the bulk critical point, and we sug- 
gest that this has repercussions for experimental studies. 

For Ising-like systems with hi = hi > in vanish- 
ing bulk field h = 0, f so iv as a function of temperature 
attains a pronounced minimum above the bulk critical 
temperature T c . The scaling function W++{y), defined 
by fsoiv/ksTc = L~ d W++(y), was determined by exact 
transfer matrix methods in d = 2 pdj |. The minimum 
occurs when y = T(Lf^) 1 ^ = 2.23, or L ~ 2.23£, 
and the amplitude at this extremum is about 6.6 times 
the Casimir value PjK ++ (0) = A ++ (d — 1). Here r = 
(T - T c )/T c and £ = ^t'" (with v = 1 in d = 2) is 
the bulk correlation length. In mean field the minimum 
occurs for L ~ 3.7£ and the amplitude is about 1.4 times 
Fy ++ (0) Borjan and Upton Q have calculated 

W++(y) in d = 3 using local functional methods and they 
find the minimum at L ~ 3.1£, and the amplitude at the 
extremum is about 2.1 times W_|_+(0). The bulk field de- 
pendence of the solvation force at T = T c is similar to 
the temperature dependence at h = 0, i.e., f so i v attains 
a pronounced minimum at some h m i n < |]l4| , ^5| . In 
d = 2 the scaling function of f so i v with h was obtained 
using the density-matrix renormalization-group (DMRG) 
method. The striking feature is that at the minimum 
this function is about 100 times bigger than the Casimir 
amplitude |l4|]! In mean field it is about 10 times the 
Casimir value [fTsf . Thus, the solvation force is much 
stronger for states which lie slightly off bulk coexistence, 
with h < 0. In this Letter we investigate the origin of 
the very attractive solvation force and argue that in a 
fluid confined by identical, strongly adsorbing walls the 
form of the near-critical, long-ranged solvation force is 
strongly influenced by capillary condensation, i.e., the 
shift of the bulk first-order transition which occurs be- 
low T c . Recall that for surface fields hi — h 2 > and 
bulk dimension d > 3, two phase coexistence occurs in 
the confined system, with finite (large) L, along a line 
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h co (T) ending in a (capillary) critical point (hcL,Tci,)- 
T CL (hi) lies below T c and h CL {hi) < [00. The 
coexistence (capillary condensation) line has a positive 
slope and is located at h < 0. f so i v exhibits a discontin- 
uous jump on crossing the coexistence line and singular 
behaviour at {Hcl,Tcl) for d > 3 the capillary 

criticality corresponds to the d— 1 universality class. The 
form of f so i v for T > Tql and h < does not seem to 
have been investigated in any detail. We perform such 
an investigation for Ising films using the DMRG method. 
Although for bulk d = 2 systems there can be no true 
phase transition for finite L, there is still a line of sharp 
(very weakly rounded) first-order transitions in the Ising 
film ending in a pseudocritical point |l^|20t|. We deter- 
mine this pseudocoexistence line and investigate f so i v for 
temperatures above the capillary condensation line and 
states approaching the bulk critical regime. 
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FIG. 1. Solvation force as a function of the bulk magnetic 
field h (both quantities in units of J) for an Ising strip of width 
L = 200, surface fields hi = h% = 10 and several reduced 
temperatures r = (T - T c )/T c ^: (+) 0.0135, (<) 0.0063, 
(V) 0.0031, (>) 0.0, (o) -0.0063, (□) -0.0126, (o) -0.0505, (A) 
-0.0546, and (*) -0.1011. 

Consider a d — 2 Ising strip with L parallel rows sub- 
ject to surface fields hi = h%{= /il) > with the Hamil- 
tonian H: 

(i) (£) 

- h/j = x a i u j + ai ~*~ ^ x X/ Ui ~*~ ^ l X] °*' 

<i,j> i i i 

(1) 

where <Ji is a spin variable taking the value ±1 and 
J > is the coupling constant. The first sum runs 
over all nearest-neighbour pairs of sites while the last 



two sums run, respectively, over the first and the Lth 
row. h is the bulk magnetic field. The total free energy 
per site of this system can be written as f(L, T, h\,h) — 
f b (T,h) + y w {T,h u h) + ±r(L,T,h u h), where f b is 
the bulk free energy, f w is the L-independent wall ten- 
sion at each wall, and f*(L) is the finite-size contribution 
to the free energy, which vanishes for L — > oo. All ener- 
gies are measured in units of J, temperature in units 
of J/ks and the width L in units of the lattice con- 
stant a. The solvation force for our system is defined 
as f so i v = -(df*/dL)x,h u h- We calculate f so i v using 
the recently developed, very accurate DMGR method, 
which provides an efficient algorithm for the construction 
of the effective transfer matrix for large systems . A 
great advantage of this method is that it works equally 
well for h ^ 0, for which no exact solutions are avail- 
able, as for h = 0. The total free energy / is obtained 
from the leading eigenvalue Ao of the effective transfer 
matrix which in the DMRG method is calculated numer- 
ically [0 . By calculating the excess free energy per unit 
area fex{L) = (f — fb) L at Lq + 2 and Lq we obtain f so iv 
by a finite difference. Figure |l| summarizes our results 
for fsoiv m a film of width L — 200 and h± = 10, con- 
ditions which pertain to the infinite surface field scaling 
limit |Q . f so i v is calculated as a function of h for several 
temperatures above and below T c . For the three lowest 
temperatures shown we find a very weakly rounded jump 
of the solvation force from zero to a negative value. As 
mentioned earlier, a discontinuous jump is characteristic 
of the solvation force at a first-order capillary condensa- 
tion phase transition [fOp^T]- At fixed large L and fixed 
temperature T < Tql , fsoiv should change abruptly from 
values appropriate to a spin down (— ): / 
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to values 

appropriate to a spin up (+): f^ olv ~ 2hm*(T), where 
m*(T) > is the bulk spontaneous magnetization [p2| . 
The calculated gradients of f so i v at the three lowest tem- 
peratures in Fig. ^ agree with the known values of 2m* (T) 
to a relative accuracy of 10 -4 . Coexistence is given by the 
Kelvin equation [^2| which for the strong surface fields hi 
considered here takes the form —h co (T) w <j(T)/Lm*(T), 
where cr(T) is the interfacial tension between coexisting 
(+) and (-) phases. The presence of thick wetting films 
of + spin in the (-) phase gives rise to non-trivial cor- 
rections which shift the condensation line to larger val- 
ues of | ft. | p3|) ; nevertheless the Kelvin equation does 
predict the correct qualitative behavior of the condensa- 
tion line at low temperatures. The Kelvin result implies 
Afsoiv = flu, - fsoiv » 2h co (T)m*(T) » -2a(T)/L, 
i.e., the magnitude of the jump should decrease in the 
same fashion as the interfacial tension, as T increases 
at fixed L. Our numerical results confirm the predicted 
trend |24| . A second signature of capillary condensation 
is a discontinuous jump (see Fig. ||) of the total adsorp- 
tion, r = *}2a = i < (J i >, from negative values character- 
istic of a (-) phase, with wetting films of + spin, to posi- 
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tive values characteristic of a (+) phase condensing in the 
slit. The simple formula Ar = T+ - T~ « 2Lm*(T) |§| 
provides a good qualitative description of the jump Ar 
obtained by DMRG for the three lowest temperatures. 
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FIG. 2. Total adsorption Fvsft (in units of J) for the same 
system and the same reduced temperatures r as in Fig. |lj 

The jumps in T are slightly less rounded than in f so iv 
but they occur at the same value of h. This is illus- 
trated further in Fig. ^ where for each T we plot the 
inflection point (*) of f so iv Note that these points lie 
on the line ho(T) defined by T = or, equivalently, by 
the maxima of the total free energy /. For tempera- 
tures higher than r = —0.0505 we observe a change in 
the behaviour of f so i v and T; see Figs, [j] and |[ As the 
temperature increases the jump of f so iv gradually trans- 
forms into a minimum, whose depth decreases and moves 
monotonically towards h — 0. r = 0.0063 corresponds to 
the temperature at which f so i v exhibits a minimum as 
a function of r for h — and L = 200 (H). At this 
temperature the minimum of f so i v as a function of h is 
still present and its depth is about 4.8 times bigger than 
I fsoiv(h = 0) |. For temperatures up to the critical 
temperature the magnetization profiles for small | h \ 
are similar to that of a "condensed" (+) (liquid) phase. 
Moreover, f so i v still varies linearly with h, as is implied 
by the approximate treatment above, although the linear 
region does shrink as r — > _ . The adsorption T exhibits 
a steep increase at r = (Fig. |^) indicating residual con- 
densation. However, the locus of ho(T) moves to larger 
| h | for t > 0, as does the inflection point of f so iv', see 
Fig. ^. The longitudinal spin-spin correlation length £y 
provides an alternative signature of capillary condensa- 
tion. If we take the transfer matrix in the infinite di- 
mension then £7 (L,T, hi,h) = — ln[Ai/Ao], where Aq 



and Ai are the largest and the second largest eigenval- 
ues. We calculated l;^ 1 as a function of h for a series 
of temperatures. At the three lowest temperatures con- 
sidered in Figs. [j] and |[ where the solvation force and 
the adsorption behave in a way characteristic of capil- 
lary condensation, ^"j -1 has a sharp minimum at some 

h m in(T) < with ^(hmin) « 0. The values of h min (T) 
lie very close to those where f so iv jumps and to ho(T); see 
Fig. ^. Recall that in Ising strips with periodic bound- 
ary conditions the two largest eigenvalues of the trans- 
fer matrix are asymptotically degenerate such that, as 
L — > oo, t;^ 1 ~ cxp(— La(T)/kBT) for states on the line 
of pseudocoexistence h = 0, T < T c . Moreover, as in the 
periodic system p5[| , our calculated S^ 1 is symmetric 
and increasing linearly in | h \ in the close neighbour- 
hood of its minimum, i.e., near pseudocoexistence. As 
the temperature increases the minimum gradually lifts 
away from zero and becomes more shallow. Its location 
remains close to ho(T) until just below T c . For T > T c , 
(T) moves towards h = 0; see Fig. H pa. 
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FIG. 3. Maxima ho(T) of the free energy (or zeros of the 
adsorption F) (solid line), inflection points of fsoiv (*) and 
minima of ^"j -1 (o) as functions of the bulk magnetic field h 
(in units of J) calculated at fixed temperature T (in units of 
J/fcs) for the same system as in Figs. [I] and ^. The bulk 
critical temperature T c ~ 2.269185 (r = 0) is denoted by 
the horizontal line. Pseudocoexistence occurs along ho(T) for 
T < 2.16. 

In the light of these various results it is natural to 
identify ho(T), at sufficiently low T, with the pseudo- 
coexistence line h co (T). Determining the pseudocritical 
temperature Tql is more difficult. However, different cri- 
teria, namely the erosion of the jumps in f so lv and in T 
and the lifting away from zero of the minimum of Ci? 1 ' 
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all yield similar estimates. We also calculated the av- 
erage "central bond energy" Ul/2 =< &l/2< t l/2+i > at 
points on the line ho(T) and dUL/2/dT on that line. The 
latter has a pronounced minimum at T w 2.145 where 
/i Q (T) w -0.00184 for the case L = 200, /ii = 10 con- 
sidered here. Although the critical point is not sharp 
in this system the various criteria do indicate a signifi- 
cant change in character for T between 2.155 and 2.160, 
providing the estimate for Tql- Such estimates are con- 
sistent with the scaling relations (ll| for {Hcl,Tcl) j p6| . 

We expect that in d — 3, where there is true coex- 
istence for T < Tcl, fsoiv and T should behave in a 
similar fashion to what is observed in d = 2, but now 
the jumps of f so iv and T will be discontinuous. Also on 
the critical isotherm T = T C l, Ar ~ (h C L - h) 1/5 , as 
h — > h~ CL , and 8 — 15, the d — 2 Ising exponent. Explicit 
mean-field results for a Landau (square-gradient) theory 
support our expectations for T < Tcl and for T > Tcl, 
where f so i v « 2hm*(T) for h — > 0~ (TsJ . Our results 
imply that the solvation force in a real confined fluid at 
temperatures near T c and at reservoir densities slightly 
below the critical value (or compositions slightly away 
from the critical composition in a binary mixture) should 
be much more attractive than the Casimir value for the 
same L, although the effect may be less pronounced than 
in d = 2. In keeping with others [^| JlO|JT2] | we suggest that 
experiments which aim to measure the Casimir force, 
e.g. by atomic force microscopy or surface force appa- 
ratus, and future theoretical work should focus on the 
r and h dependence of f so iv, i.e., the scaling function, 
not just the Casimir amplitude. Finally we note that the 
existence of a long-ranged, strongly attractive solvation 
force between two large colloidal particles immersed in a 
near-critical fluid can have ramifications for aggregation 
or flocculation of the particles |l2| . Studies such as the 
present should indicate where aggregation is potentially 
the strongest. 
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